LIMITS ON THE DOMAIN OF COUPLING CONSTANTS 
FOR BINDING iV-BODY SYSTEMS 
WITH NO BOUND SUBSYSTEMS * t 



Jean-Marc Richard*- 1 -*''- 2 ) and Sonia Fleck^) 
^ Institut des Sciences Nucleaires 
53, avenue des Martyrs, 38026 Grenoble, France 
( 2 ) European Centre for Theoretical Studies in Nuclear Physics and Related Areas (ECT*) 
Villa Tambosi, via delle Tabarelle 286, 1-38050 Villazzano (Trento), Italy 
( 3 ) Institut de Physique Nucleaire 
4-3, boulevard du 11 Novembre 1918, 69622 Villeurbanne, France 

(February 9, 2008) 



Abstract 

We study the domain of coupling constants for which a 3-body or 4-body 
system is bound while none of its subsystems is bound. Limits on the size of 
the domain are derived from a variant of the Hall-Post inequalities which re- 
late iV-body to (N — l)-body energies at given coupling. Possible applications 
to halo nuclei and hypernuclei are briefly outlined. 
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For the sake of clarity, we shall define in this Letter a "halo" as a 3-body quantum 
system that is bound while none of its 2-body subsystems has a discrete spectrum. More 
generally an iV-body halo is bound while none of its subsets is stable against spontaneous 
dissociation. This is more restrictive than the usual meaning of a weakly bound system with 
a very extended wave function. These systems are sometimes called "Borromean" after 
Borromean rings, which are interlaced in such a subtle topological way that if any of them 
is removed, the other two would be unlocked. 

Halo states are seen in nuclear physics fl|||. For instance the (a,n,n) system is bound 
( 6 He), while (a, n) and (n, n) systems are both unbound. There is a cooperative effort of all 
attractive potentials to achieve the binding of 6 He. 

The halo phenomenon shows up in simple potential models, as seen from explicit cal- 
culations on specific isotopes or from the rich literature on the related Efimov effect 
H or Thomas collapse ||, which prove that a 3-body system is, indeed, more easily bound 
than a 2-body one. 

Consider a short-range potential gV(r) acting between two particles of mass m separated 
by r. Even if V(r) is attractive, a minimal strength g^/m is needed to achieve binding, where 
#2 is independent of m. For instance, in a Yukawa potential V = — (exp(— /xr))/r, one can 
fix the energy and distance scales so that m — \i — 1 without loss of generality, and one finds 
#2 — 1.680 ||. (A simple argument by Dyson and Lenard 0] shows that g% > v2.) Now, if 
one considers three identical bosons with mass m — 1 interacting through XX ex P — r ij)/ r iji 
where is the relative distance between particles i and j, one can look at bound states by 
variational methods or by solving the Faddeev equations, and one finds binding for g > g%, 
with g\j g% — 0.804. We call this a 20% window for halo binding. 

If one repeats the above calculations for other short-range potentials, one finds similar 
values for g^/g^, for instance 0.801 for an exponential form V(r) = — exp(— r), or 0.794 
for a Gaussian V = — exp(— r 2 ). Such quasi-universality, already noticed in §|, is not too 
surprising. In the weak coupling limit, the wave function extends very far away and does not 
really probe the detailed structure of the potential. Any short-range interaction is almost 
equivalent to a delta function in this limit. 

It seems therefore very likely that the window for halo phenomena is limited, i.e., g%jg\ 
cannot be made arbitrarily small by tuning the shape of V(r). We shall show below that 

£>- a) 

for any potential, and we outline how to derive similar inequalities for asymmetric 3-body 
systems, or for 4-body or more complicated systems. 

The method of deriving ([!]) is directly inspired by the Hall-Post inequalities |§, which 
have been recently rediscovered |)lC)f , applied to hadron spectroscopy [JO]] an d generalized to 
the case of unequal masses |12| . Earlier applications were mostly motivated by considerations 
on the stability or instability of matter. The Hamiltonian for three bosons of mass m can 
be written as 

i=i z '' n i<j 

6m ^3m\ 2 J y K l ° h V ; 
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i.e., introducing the translation-invariant part of the Hamiltonian H 



N 



2 (3m/2,^) = -2^ 
i<j i<j 



iJ 3 (m,s)=£#? } (3m/2 



^E^'K3 9 /2). 



(3) 



Hence, from a simple application of the variational principle, the ground-state energies 
for a given (large enough) coupling g satisfy 



E 3 (m,g) > 3E 2 (3m/2,g) = 2E 2 {m,3g/2), 



(4) 



which is the simplest form of the Hall-Post inequalities. The decomposition (f|) also implies 
that if if 3 has to support a bound state, each H 2 should produce a negative expectation 
value in the corresponding wave function, and thus 3g/2 > g 2 , q.e.d. 
A straightforward generalization of (|f) to iV bosons is [10 



H N (m,g) 



N 



N rw 



k=l 



f Nm 
\N - 1 



(5) 



where the superscript in H means that the fc-th particle is omitted. Saturating with the 
ground state oi H N gives 



9n> 



N-l 



N 



~9n-ii 



(6) 



i.e. Ng c N increases with N. 

For numerical illustration in the N = 4 case, we adopted a variational method that is 
widely used in quantum chemistry fl3| . It is based on trial wave functions of the type 



E 



exp 
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X; 



(7) 



'•j 



where {xj} is a set of relative Jacobi coordinates. Symmetry is properly implemented by 
imposing relations between neighbouring coefficients c^ n ' and definite-positive matrices af< n \ 
After numerical minimization, we obtain g c J g 2 — 0.67 for a Yukawa potential, i.e. a 13% 
window for a genuine 4-boson halo, once gl/g 2 is subtracted from g\j ' g 2 . 

Note that the bounds (p and @ are not expected to be saturated, since the Hall- 
Post inequalities become equalities only for harmonic oscillators, which are far from the 
short-range potentials we consider here. 

Similar inequalities can be written down in a variety of situations. Let us give some 
examples. Consider first N identical particles with mass m = 1 in the field of an infinitely 
massive source. The Hamiltonian is defined as 



H 



N 



N 

E 



2m 



+ hV{ri 



+ <?E^(' 



(8) 



The short-range attractive potentials V and W can be normalized such that a single particle 
is trapped around the source for h > 1, and two particles are bound together for g > 1. The 
non-trivial domain is thus to be found inside the square (h < 1, g < 1). 
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Let h,N{g) be the boundary for binding N particles around the source. This means that 
halo type of binding occurs between h^^g) and fiN-i(g)- One expects all curves to merge at 
A(h = 1, g = 0) as g — > 0, since the particles become independent in this limit. If one takes 
the expectation value of the identity 



1 



H 3 (h,g) = -J2H^\h,2g) 

i<j 



within the ground state of H 3 , one gets 



h 3 (g) > h 2 (2g), 



and more generally 



Ms) > h N-i {jj~r^9 

For N = 2, one can write the simple decomposition 

I - 



Ho 



a 



p{ + hV( n ) 



a 



p{ + hV{r 2 ) 



a 



pi + P2) + gw{r 12 ) 



(9) 



(10) 



(11) 



;i2) 



for any < a < 1. To get (H 2 ) < 0, one needs at least one of the square brackets having a 
negative expectation value. This excludes the triangle {h < a, g < (1 — a)} shown in Fig. |l|. 

Two remarks on this simple lower bound are in order. First, the actual boundary is 
expected to be concave. The couplings h and g enter the Hamiltonian linearly, so if the 
minimum E of H 2 vanishes at both P(h,g) and P', one has [ 14] 



E(XP + (1 - X)P') > \E{P) + (1 - X)E(P') = 0, 



(13) 



for any < A < 1. 

Secondly, while the limit A(h — 1, g — 0) of truely independent particles obviously 
belongs to the boundary, B(/z = 0, g — 1) is likely to be in the continuum, since a weakly 
bound (1,2) system needs a minimal attraction h to remain trapped by the source. 

A more elaborate decomposition leads to an improved boundary which better complies 
with the above remarks. We provisionally restore a finite mass M for the third particle and 
write as in 
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H { M) 



222 
^ + ? + 2M + kV{rn) + kV{r23) + gW{ri2) 



Pi + P2 + P3, 

1 + X 



+ E 

i=l,2 

+ 



[Opi 

1 + 2x7 \2 + M 



bp2 + b'p 3 ) 



[i + 2 X y 



x(l + x) — 



2M 



Pi ~ ZP3 

1+X 



Pi ~ P2 



+ hV(r l3 ) 



+ gW(n 2 ) 



(14) 



where the momentum (p, — xps)/(l + x) is the conjugate of the relative distance r i3 , and b 
and b' are known functions of M and x. In the limit M — > 00, we read off from (fLU) that 
H 2 would never support a bound state if 
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ft, < 1 and g)— J — < 1. (15) 

\l + xj Ax(l + x) 

This is the inner part of the parabola shown in Fig. [I], from which we get a crude lower limit 
on the minimal coupling h 3 (g) to bind three bosons around the source, as per Eq. (|T0[) . 

The decomposition (|l4]) can be used for finite M. Consider for instance the case where 
M = 1 and g = 0. One should restrict to x(l + x) > 1/2 in order not to introduce a negative 
reduced mass. The two particles, which do not interact with each other, can be bound 
simultaneously below the critical coupling h = 1 for individual coupling. Each particle 
benefits from the increase of the reduced mass provided by the other. However, by choosing 
the optimal parameter x in (0), one can easily deduce that the window for halo is limited 
to h > 1/2 + y/3/4 ~ 0.93, i.e. at most 7%. 

In a situation where the three masses or the three couplings are different, the most 
general decomposition involves two parameters instead of the single x in (14). So the 
analysis becomes slightly more involved. 

We now write the Hamiltonian for two identical particles of mass m, and two others of 
mass M 

p 2 P 2 

H ± = ^ + Y ^ + #12^12 (?"12) +#34 ^34(^34) + gmAlYl^Mirij), (16) 

i=l,2 i=3,4 i,j 

including up to three different potentials. It can be rewritten as 
Hi = (pi + p 2 + P3 + P4) • (6pi + bp 2 + b'p-i + fe'p 4 ) 

+«i2 ( Pl 2 + 9i2Vi 2 (r 12 ) + a u ( P3 2 P J + 93iV 34 (r u ) 

+ «[«Pi - (1 - a )Pj] 2 + gmMV mM {rij) (17) 

For given a, one can solve for b and b', as well as for the inverse masses, with the result 

Ol2 = — (1 - « 2 ) - TTa 2 , 

m M 

a u = --(l-a) 2 + -J-a(2-a), (18) 

1 1 

+ 



4m 4M 

If one rescales the couplings to the critical value for 2-body binding with the appropriate 
inverse reduced mass, m~ l for gi 2 , M~ l for g^, and (m~ l + M~ 1 )/2 for g m M, then cannot 
support a bound state if simultaneoulsy 

gn < 1 - a 2 - (m/M)a 2 , 

934 < -(M/m)(l - a) 2 + a{2-a), (19) 

g m M < 1/2. 

Interestingly, the condition on g m M decouples. Hence the domain for 4-body binding without 
2-body binding consists at most of 1/2 < g mM < 1, and, in the (gi 2 , g 34 ) plane, the area 



5 



between the unit square and a parabola, as shown in Fig. ||. One should of course exclude 
the domain corresponding to 3-body binding to get a genuine halo. 

Our interest in 3-body systems was clearly triggered by nuclei like 6 He or n Li with two 
neutrons weakly attached to a compact nucleus [Q]. We assume a spin singlet state for the 
two neutrons, so that their spatial wave function is symmetric. States with more than two 
neutrons in the halo unfortunately escape direct application of our results, because of the 
Pauli principle. Our result on (m, m, M, M) configurations is perhaps relevant for (n, n, A, A) 



hypernuclei with strangeness S = —2, a field of intense theoretical studies jlSj . 
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Brac, and T.T. Wu. The Institut de Physique Nucleaire is supported by Universite Claude 
Bernard, the Institut des Sciences Nucleaires by Universite Joseph Fourier, and both by 
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FIGURES 




FIG. 1. Lower limit for the domain of halo binding of two bosons in the field of a static source. 
The coupling to the source is denoted h, while g is the interparticle coupling. Two-body binding 
occurs for h > 1 or g > 1 . The straight line AB is deduced from the simple decomposition (|l2|) , the 



parabola from (14). The dotted parabola is a lower limit for binding three bosons in this central 
field, as deduced from 
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FIG. 2. Lower limit for binding four particles with masses (m, m, M, M). The couplings g\2 and 
534 are normalized to the critical coupling for binding (m, m) and (M, M), respectively. Meanwhile 
the (m, M) coupling should be at least half the critical coupling for binding m to M. A value 
M/m = 2 is assumed here for the drawing. 
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